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1. Introduction 

Symmetry is an equivalence of different physical situations [ ]. Such an equivalence 
manifested in an invariance of a certain set of observables can be formalized in terms of 
commutation relations between observables and Hamiltonian of the system in question. 
The existence of conserved quantities, also having no classical counterpart, extends 
an amount of an information accessible for a given quantum system. In an extreme 
cases one can meet analytically solvable models such as harmonic oscillator, Jaynes- 
Cummings model or hydrogen atom. The more symmetries are recognized (together 
with related conserved quantities) the more knowledge about the system is accessible 
for the researcher. It is a kind of truism in studies of quantum systems. 

In the paper we consider a quantum model consisting of a two-level system (qubit) 
interacting with a single mode bosonic field (representing e.g. electromagnetic radiation) 
with frequency u. Hamiltonian for that system is assumed to be 

H = (3a:, + Acr^. + ua'a + cr^ (g) {g*a + go)) , (1) 

where a and are the creation and annihilation operators of the bosonic field i.e., 
[a, a^] = I, whereas cr^ and denote the two of the Pauli spin matrices. The term 
Paz is an unperturbed energy of the qubit, with possible eigenenergies ±/3, while 
Acr^ describes a tunneling of the two-level system in the absence the bosonic field 
(spontaneous transition). Finally, the coupling constant g reflects the strength of the 
interaction between the systems. 

The above Hamiltonian is the well known Rabi model [2, 3]-probably the most 
influential model describing fully quantized interaction between matter and light. 
Although, the model originates from quantum optics [ ], its applications range form 
molecular physics [ ], solid state (see Refs. in [{i]), to the recent experiments involving 
cavity and circuit QED [7, S]. This model can be implemented by means of rich variety of 
different setups such as Josephson junctions [ ], trapped ions [UJ], superconductors [11] 
or semiconductors [ ], to name a few. 

Despite simplicity of the model its Hamiltonian (1) cannot be diagonalized exactly 
when A 7^ 0. Although, some progress has been reported recently [13, 14], the exact 
analytical formulas for eigenvalues and eigenfunction are still missing. There is a wide 
spectrum of available approximation techniques, including RWA [i i] (leading to the 
famous Jaynes-Cummings model [' ]), which allow approach the problem from many 
different directions. 

At this point a natural question concerning the existence of any symmetries in this 
model (and related constants of motion) arises. The presence of a constant of motion 
within an unsolvable model is not only of great theoretical importance, but also can be 
applied to build its experimental realization [17]. If /3 = the Hamiltonian (1) remains 
unchanged for a, — J- — o"z and a ^ —a (hence — )■ —a^). The symmetry operator Jq, 
which generates this transformation {e.g. fulfills [H, Jq] = 0) simply reads Jq = ® P, 
with P = exp(i7ra^a) being the bosonic parity [ ]. This is the well known result: Still 
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unsolvable system possesses a symmetry when /3 = 0. 

On the other hand, when /3 7^ one can leave H unaffected by simultaneously 
changing cr^ — )■ — cr^, a — > —a and in addition /3 — )■ —f3. This instantly brings a question: 
How does the corresponding symmetry generator J looks like? Unfortunately, this 
question has not been answered so far. Moreover, quite recently it was conjectured [13] 
that the Rabi model (1) does not posses any symmetry at all (except the trivial one, 
related to the total energy conservation) as long as /? 7^ 0. In other words, it was 
expected in [l .>] that the only self-adjoint operator J such that [H, J] = is J = H. 

On the basis of the results reported here, we prove that this conjecture is not true. 
We show that in the Rabi model one can find a self-adjoint involution J, that is = Ib, 
such that HJ = JH. We also offer arguments why this constant of motion is so hard to 
find in practice, in comparison with the case of unbroken Jo symmetry. Possibilities of 
the exact diagonalization of (1) are also discussed. 

2. Main result 

We begin with formal rewriting the Rabi Hamiltonian (1) as a matrix with operator 
entries: 

where H-t := ua^a ± (^g*a + ga^^ ± /3. (2) 

Customary, the parameters A and (3 denote AIb and /3Ib, respectively and Ib stands 
for the identity on the bosonic Hilbert space Hb- 

The matrix representation (2) of the Rabi model is established via natural 
isomorphism ® "Hb ~ "Hb © ^b- Such an identification is usually invoked in order 
to simplify purely algebraic calculations [19]. This is not the reason why we use this 
idea. Instead, we are going to attack the problem in question by using concepts of block 
operator matrix [ ] along with its relation to the operator Riccati equation [21-24]. 

Before, we need to clarify some technical aspects concerning the matrix (2), such 
as e.g., its domain I'(H), so that it may become meaningful mathematical object. 
Unfortunately, in many papers addressing physical aspects of the Rabi model this is 
not a primary issue. It is needless to say, that one cannot then take advantage of very 
powerful existential mathematical theorems (like e.g., the famous Banach fixed point 
theorem [2")]), simply because is not known whether the premises of these statements 
are met. 

As a first step toward constructing J, we define domains := V{B.±) on which 
both operators }l± are self-adjoint. Since the off-diagonal elements of H are bounded, 
we conclude that also H* = H on 'D(H.) = T>+ © !)_. To end this, a proper definition 
of creation and annihilation operators is required. Both a and are unbounded [26] 
and hence the canonical commutation relation can only holds on some (dense) subspace 
T>2 of l-L-Q- Let us assume that Vi is a dense set on which both a and are adjoint, 
that is {0))* = a and a* = a"^ . At this stage, the existence of Vi and V2 having the 
desire properties is not obvious. Result of detailed construction given in [ ] can be 
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summarized as follows 

{oo oo ^ 

Y,Un)enB: ^n'=|en|'<ooL k = l,2, (3) 

n=0 n=0 J 

where {|^^)}J^o canonical (orthonormal) basis in /2(~ "Hb)- This is rather obvious 

choice since a, and a^a need to produce normalizable states. Having (3) in place, we 
define (see also [28, 29]) 

oo oo 

a\ij) ■.= Y,V^Un-l), a^i/j) := ^ V^i + l^nln + 1), |^) G Pi. (4) 

n=l n=0 

From (4) it immediately follows that a^|n) = i/n + l|n + 1) and a|n) = — 1). In 

most textbooks on quantum mechanics the latter relations serve as the very definition 
of the creation and destruction operators. Technically, this is not a good choice for a 
definition due to the problem with closeness. This can be overcome by taking closures 
given by (4) instead. 

A very basic result from operator theory (see e.g., Theorem 4.2.7 in [.>()]) states 
that if A is closed on T'(A) then A*A is a positive self-adjoint operator. Moreover 
'D(A*A) C 'P(A) is a core of A, that is A is a closure of its restriction A|x>(a*a)- Since 
the operators Hj- can also be written as 

H± = a;fa±^V fai^) ±/3-^, on V2, (5) 

they are both self-adjoint and their common domain is a core of both a and 
. Therefore, the total Hamiltonian (1) is well defined self-adjoint operator on 

After resolving technical nuances of the Rabi model we introduce quadratic second 
order operator equation: 

AX^ + XH+ - H_X - A = 0. (6) 

Within the framework of the operator matrices theory, Eq. (6) is known as the Riccati 
equation associated with matrix (2). As we will see, most of the relevant problems 
regarding the Rabi model (1), including possibility of its exact diagonalization, can be 
reduced to the mathematical questions concerning solvability of this equation. 

Because (6) is an operator equation its solution can be understood in either a week 
or strong sense. There is also a third possibility-an operator solutions [ ], which is not 
discussed in this paper. Let us briefiy clarify this notions. 

Definition 1. A hounded operator Xq on H is said to be a weak solution to the Riccati 
equation (6) if 

A{Xl(j),i;) + {XoH^(P,ij)-{Xo(P,H,ij)-A{(f^,i;) = 0, for ,\4>) e V^. (7) 



New symmetry in the Rabi model 



5 



A bounded operator Xq on "H such that Ran{XQ\'p^) C 1^2 one? 

/\Xl\^) + XoH^\,p)-H_X,\,p)-/\\,p) = Q, for \^)eV^, (8) 

is called a strong solution of (6). 

Of course, a strong solution is also a week solution. From technical point of view 
it is often easier to prove the existence of a weak rather than strong solution. However, 
in quantum mechanics we are usually seeking for strong solutions of operator equations 
as this is how we are used to equate the operators. Fortunately, this two notions are in 
fact equivalent [32]. Nevertheless (6) there is no general method allowing to find either 
its weak or strong solution. For this reason, the following theorem offering criteria of 
solvability is of great importance for us. 

Lemma 1. Let H± be possibly unbounded self-adjoint operators on domains V{H±) in 
separable Hilbert space Ti. Also assume Vi ^ and V2 are bounded operators onl-i. If 
the spectra (j{H±) are disjoint, i.e., 

d := dist{a{H+),a{H_)) > 0, (9) 

and if one imposes the smallness assumption on operators Vi, V2 of the form 

V\\yi\\\\v2\\ < -, (10) 

TT 

then the Riccati equation 

XViX + XH+ - H_X - V2 = 0, (11) 
has a unique weak solution Xq in the ball 

X e B{n) : \\X\\ < (12) 

The weak solution satisfies the estimate 



l^oll < ^ I ll^illll^2|| I . (13) 



In particular, if 



2d 

l^i|| + ||^2||<-, (14) 



TT 

then Xq is a strict contraction, that is, \\Xq\\ < 1. 

An elegant and compact proof of this statement, based on the Banach fixed point 
theorem, can be found in [ ]. Now we can prove our main result. First we show that 
if exists a solution to (6) then there is also an operator generating symmetry in the 
system (2). Second, we show that under certain conditions imposed on the parameters 
a, /3, u Eq. (6) is weakly solvable. 
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Theorem 1. Assume there exists a weak (hence strong) solution Xq to the Riccati 
equation (6). Then, there is also a self-adjoint involution J such that JH = HJ, where 
H is given by (2). Moreover, in terms of Xq the generator of the symmetry J reads 



XqJqX*q - 



where Jq = 2(1 + X*qXo) 



Proof. Let ^(Xq) be the graph of Xq, that is 



6^(Xo 



Xn 



(15) 



(16) 



Since Xq is a strong solution to (6), XqIV^) G V2 and Xo(H+ + AXq) IV") = (H_Xo + A)|^) 
for 1-0) G 1^2 • Therefore 
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' m ' 
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H 








nV2 
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(H+ + AXo)|' 
Xo (H+ + AXo) 



Ge?(Xc 



(17) 



^(Xo)"*", which in this case reads 
^(Xo)^ = 



-X* 



G ?/b ©^B : IV^) G Hi 



is H-invariant as well. Since Xq is bounded its graph forms a closed subspace of 7^b©'Hb 
and hence the decomposition "Hb © "Hb = ^(Xq) © ^(Xq)"*" holds. Therefore, each state 
|\E') G ^^(H) of the composite system can be uniquely decomposed as |\E') = |\E'i) © {^2)) 
where |^i) G ^(Xq) and (*2|^i) = 0. 

Let P+ be the projection onto ^(Xq), then P+H|^i) = H|^i) and P+Hj^s) = 0, 
thus assuming, for a moment, that P+X'2 C 1^2 we obtain 

H(P+|*i)©P+|^2)) =H|^i), and P+ (H|*i) © H|^2)) = H|^i). (19) 

Therefore HP+|^) = P+H|^), for all |$) G V2. The inverse (1 + X*Xo)-i exists and it 
is a bounded self-adjoint operator on Hb and therefore P+ can be expressed as 



JoXq 



(20) 



p^^i Jo 

2 XqJo XqJoXq 

It is a matter of straightforward calculation to see that (20) indeed projects onto ^(Xq). 

Since J = 2P+ — 1 the only question we need to addressed in order to conclude 
the proof is whether P+|^E') is again in I^(H) for G I^(H). Because Xq is weak (and 
hence strong) solution to Eq. (6) we have X0P2 C T'2 and moreover, for every |0) G T>2, 
the function /('0) := (H+'0,Xg0) is continuous on T'2- Indeed, from (7) it follows 

|2 



1/(^)1 <M^||V^||, where = a||0|| UXof + ||H_0|| ||Xo|| + a||0||. (21) 

Since / is continuous as just proved we obtain Xq|0) G ©(I!*) = V2, thus also 
XqT'2 C T'2- As a result ^q^T>2 C T'2 and because Jg ^ is invertible this in fact means 
JoT'2 = T'2. The latter implies that P_,.T'(H) C T'(H) which concludes the proof. □ 
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Theorem 2. Let f3, u, A and assume the following conditions 

and ->-. (22) 
uj P 2 

Then there exists a unique weak (hence strong) solution of Eq. (6) such that \\Xo\\ < 1. 
Under assumptions (22) the Rabi model (1) has a self-adjoint symmetry given by Eq. (1). 

Proof. Since V^; = i(xa^ — x*a) is self-adjoint for x G C, the unitary Weyl operator 
Dj; = exp(iVa;) is weU defined. Moreover D* = D.^^ and therefore we have the well 
known result: 



H± = (^coa^a ±(3- D, 
By virtue of a'^aln) = n\n), for n G N we obtain 



(23) 



(j(H±) = [un ±/?-^:nGNl=a;NU {±/3} - (24) 

[ U ) CO 

We see that if 2/? is not multiple of u then the distance 

dist(a(H+), a(H_)) = inf{|w(n - m) + 2(3\ : n,m E N} = 2f3. (25) 

and therefore the spectra o'(B.±) are disjoint, i.e. (9) holds. In addition, both the 
smallness assumption (10) and (14) imposed on the off-diagonal elements are satisfied 
as long as 2A > tt/S. Therefore, according to lemma (1) there is exactly one solution Xq 
to the Eq. (6) and ||Xo|| < 1. 

The second statement of the theorem follows immediately form theorem (1). □ 



3. Discussion 



First let us take a closer look at the /3 = case. In this situation the spectra (24) 
overlap, hence the separability condition (9) is not satisfied. Therefore, one cannot 
invoke theorem (1) to establish the existence of a solution to (6). However, in this 
particular case the spectra cr(II±) are identical and }i± can be transformed one into 
another by the same bosonic parity operator which generates the symmetry Jq. As we 
shortly see this is not an accidental coincidence as P happens to be a solution to (6). 
Indeed, first note that 

oo 

P|^) = E(-1)"^«I^)' ^ith en = (n|^). (26) 

n=0 

This equality can serve as the definition of P. 

From (26) it immediately follows that P is bounded and in addition to that also 
Ran(P|x)2) C ©2 (if n^n are square-summable so are {—l)"'n$,n), which together with 
Eqs. (4) gives PaP = —a as well as Pa^P = —a^ and finally it leads to PH-|-P = B.^. 
Since P is also a self-adjoint involution, it solves (6) as stated. This solution is not 
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unique; indeed, — P satisfy Eq. (6) as well. In this case the symmetry generator (15) 
form theorem (1) is equal to Jq (or —Jo) as one should expect. 

If condition (22) are met, in particular for /3 7^ 0, the spectra H± are separated 
and Eq. (6) possesses exactly one solution which according to (1) corresponds to a 
symmetry generator J of the Rabi model. The only problem is that this solution is 
unknown. However, one can attempt to simphfy the problem by putting Xq = Y/jP; 
then 

aY^PY^ + [Y^, + 2/3Y^ - aP = 0, (27) 

where we redefined H+ so that it reads (5) with = 0. When (3 = this equation 
becomes trivial and the solution can be retrieved to be Yq = 1. On the other hand, 
as long as /3 = 0, under (22), the Eq. (27) satisfies premises of the lemma (1), hence a 
unique Y^ exists and ||Y^|| < 1. Note, if the inverse Y^^ exists as well then 

aY/PY^i + [Y^\ H+] + 2(-/3)Y^i - aP = 0, (28) 

and therefore Y_^ = Yj^. Although, we cannot solve this equation either, the latter 
equality indicates to what class Y^ belongs. One can also verify that operator Y^ is 
neither self-adjoint nor anti-self-adjoint (Y^ = — Y^). Indeed, assuming to the contrary 
that Y^ = Y^ exists, it then would imply the following separation into a self-adjoint 
and anti-self-adjoin part: 

aY/^FYp + 2/3Yi3 - aP = 0, and [Y^, H+] = 0. (29) 

Both this equations can be solved separately, but the solutions does not agree which each 
other unless (3 = 0. Similar arguments shows that the condition Y^ = — Y/3 necessary 
leads to Y^ = what contradicts (27) even when (3 = 0. 

The solution to Eq. (6) not only gives rise to the symmetry generator in the 
Rabi model, but also can, in principle, be explored to obtain its eigenfunctions and 
corresponding eigenvalues. Let us briefly discuss this concept. 

Since both ^(Xq) and ^(Xq)^ are H-invariant, thus if is an energy eigenstate 
it follows that either G ^(Xq) or G ^(Xq)"*". We can say even more than that. 
Namely, let us define, on V2, closed operators: Z+ = H+ + AXq and Z_ = H_ — AXq, 
then form (17) it follows that for all eigenstates of H belonging to ^(Xq) 



where Z+|^a) = A|^a). (30) 



Using similar arguments one can show that all eigenstates from ^(Xq)"*" are of the form 



-xSI^a) 

10a) 



, where Z_|0a) = A|0a). (31) 



It can be proven that Z± are self-adjoint on Hilbert spaces {(Hb, ((1 + XqXq)-, •)) and 
(?^B,((1 + XoXo)-,-)), respectively [20]. Moreover, f7(H) = a(Z+) U (t(Z_) and the 
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following similarity relation holds 
extending the notion of parity chains introduced in [13]. 





a" 


S = 


H+ + AXo 


A 


H_ 




H 



where S 



1 -XI 

Xn 1 



(32) 



4. Summary 

We have recognized a novel symmetry in the Rabi model and constructed its generator 
J. Although, this symmetry is nonlocal, unlike e.g. Jq = ® e"^"^^", it is a self-adjoint 
involution and therefore can be considered as a generalized parity of the Rabi model. 
Invoking physical nomenclature, the Rabi model is invariant with respect to this parity 
or it has an unbroken generalized Z2 symmetry. The latter terminology is often used 
in a different (local) context when it is stated that /3 7^ corresponds to the broken Z2 
symmetry as in this case [H, Jo] 7^ 0. Our aim was to generalize the notion of local 
parity constructed for total system (combined by parity operators of the individual 
subsystems: and e™"' ") to the nonlocal one for /3 7^ 0. 

Our results are not of purely existential character. By means of a solution to 
a quadratic operator equation-the Riccati (operator) equation, we have derived an 
explicit formula for the symmetry operator J and formulated conditions (22) (range 
of parameters) guaranteeing its existence. The question whether symmetries in the 
Rabi model can exists under conditions other that (22) is a subject of current intensive 
investigations. 

Finding solutions to the Riccati equation can be related to the question concerning 
solvability (diagonalization) of the Rabi model. This equation offers a 'measure of 
difficulty' of the original problem. Although, our analysis was mainly focused on the 
Rabi model, the presented technique can in principle be extended to general two-level 
open quantum systems. 
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